Directional Emission from an Optical Microdisk Resonator with a Point Scatterer 
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We present a new design of dielectric microcavities supporting modes with large quality factors 
and highly directional light emission. The key idea is to place a point scatterer inside a dielectric 
circular microdisk. We show that, depending on the position and strength of the scatterer, this 
leads to strongly directional modes in various frequency regions while preserving the high Q-factors 
reminiscent of the whispering gallery modes of the microdisk without scatterer. The design is very 
appeaUng due to its simplicity, promising a cleaner experimental realisation than previously studied 
microcavity designs on the one hand and analytic tractability based on Green's function techniques 
and self-adjoint extension theory on the other. 

PACS numbers: 42.55.Sa, 42.25.-p, 42.60.Da, 05.45.Mt 



Introduction. — Modern fabrication techniques allow 
one to build dielectric optical resonators on a microscopic 
scale. Light is trapped by utilizing the principle of total 
internal reflection. These microcavities have great poten- 
tial for a wide range of applications and studies in laser 
physics and microphotonics pj, , like the realization of 
miniature laser sources, the creation of dynamical filters 
for optical communications and the suppression and en- 
hancement of spontaneous emission. 

For many practical purposes it is crucial to have mi- 
crocavities that possess resonances with long lifetimes 
(which is a prerequisite for low threshold lasing) and 
highly directional emission patterns. The lifetimes are 
characterized by the so called Q-factor given by Q = 
cj/Acj, where uo — iAa;/2 is the complex frequency of 
the resonance with Ao; being the linewidth or inverse 
lifetime. The best known example of modes with high 
Q-factors are the so called whispering gallery modes 
(WGMs) which are the optical analogues of the acous- 
tic waves evolving along the walls of convex shaped halls 
first studied by Lord Rayleigh in the 19th century. The 
experimental realization of a thin microdisk laser based 
on WGMs was first reported in [3]. Theoretical studies 
0] show that the WGMs of an ideal circular microcavity 
can lead to very high Q-factors of the order 10^ — 10^^. 
Despite a significant degradation due to imperfections on 
the disk boundary, inhomogeneity of the refractive index 
inside the disk, effects of coupling to the substrate etc, 
typical experimental Q-factors of the dielectric disk reso- 
nances remain quite high, usually ^ 10^, with the record 
value of Q 5 X 10^ ^5[. However, the applicability 
of circular microdisk lasing cavities is limited by their 
isotropic light emission. In order to obtain a directional 
optical output one has to break the rotational symmetry, 
for example, by deforming the boundary of the cavity 
@, 0, Q . This significantly improves the emission direc- 
tionality but typically spoils the Q-factors. Another ap- 
proach to breaking the symmetry is to insert an obstacle 
like a linear defect [9| or a hole [lO| into the microdisk. 
This indeed allows one to obtain resonances with very 



large Q-factors and relatively high directionality. How- 
ever, all the systems mentioned above require quite ex- 



tensive numerics like the boundary element method [11 ] 
or the ^-matrix approach 0, |l2[ to find resonances and 
optimize design parameters. 

In this Letter, we propose a much simpler method 
which significantly improves directionality of the modes 
of conventional microdisk resonators while keeping their 
Q-factors high (> 10^). The symmetry is broken by plac- 
ing a point scatterer within the inner region of the mi- 
crodisk, see Fig. [H It turns out that such a geometry 
improves the emission directionality of microdisk modes 
for a wide range of frequencies, especially in the visible 
spectrum. Moreover, this approach is to a large extent 
analytically tractable enabling a systematic optimization 
of the design parameters (location and strength of the 
scatterer) with only modest numerical effort. Closed sys- 




FIG. 1: Microcavity of refractive index n and radius R with a 
point scatterer at the distance d from the center. The external 
medium has refractive index riext = 1. 

tems with a point scatterer have been extensively studied 
in the context of 'quantum chaos' [13] showing that the 
spectral properties of such systems can be obtained from 
self-adjoint extension theory [14]. For an open microdisk 
we will mainly refer to these results. Similarly to the 
case of closed systems, it turns out that the resonance 
wavefunctions of the open microdisk with a point scat- 
terer are essentially given by the Green's function of the 
microdisk without scatterer (unperturbed microdisk). 

Unperturbed Microdisks. — If we treat a microcavity 
as a passive object, we can find its resonances from 
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Maxwell's equations with a refractive index independent 
of the EM field. For a zero axial momentum EM field, i.e. 
for waves with /c^ = 0, where z is perpendicular to the 
disk plane {xy\ a thin passive microdisk can be mod- 
eled as a 2D dielectric disk of radius i?, with the effective 
refractive index neff(r) = n, which takes into account 
the material as well as the thickness of the microdisk. 
In this model, Maxwell's equations reduce to two scalar 
Helmholtz equations corresponding to TM and TE po- 
larizations, respectively. 

In this Letter we consider only TM modes. The elec- 
tric field is then of the form E = Ez{x^y) and for a 
wavenumber k = uo/c^ Ez satisfies the Helmholtz equa- 
tion (V^ k'^n^)Ez = 0. Due to circular symmetry 
this equation can be separated in polar coordinates (r, (j)) 
where the TM modes are then characterized by azimuthal 
and radial quantum numbers m = 0, ±1,±2,... and 



q = 1,2,3,..., respectively. The corresponding Green's 
function is given by (see, e.g., [lH) 

G(r,ro,fc)= Yl 2nr W E^(r<'k)E2{r>,k) , 

ra— — oo 

where £^1^2 (^) are solutions of the homogeneous equation 

W = W{Ei^E2) is the Wronskian evaluated at ro, and 
r< (r>) is the smaller (larger) of r and ro- The physical 
boundary conditions require Eiir) to be finite at r = 
and E2{r) to be an outgoing wave for r ^ 00. Requiring 
£^1^2 to be smooth dX r = R then leads to 



G(r,ro, A:) = < 



jHQ(kn\Y -yq\) ^Crn cos [m {cp - (po)]Jrn{knr<:) Jrniknr^), r<,r^ < R, 



A, 

m=0 
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— — ^ -—ernCos[m{^-^Q)]Jrn{knr<)Hrn{kr^), r< < i?, r> > R, 



(1) 



2TikR A 

m=0 



where = 2 if m 7^ and = 1 if m = 0, and 

Ajn = nRrn {kR) Jm (kuR) — {kR) Jra (kuR) , 
Cm = Hm {knR) {kR) - nH^ {knR) {kR) . 

The functions Jm and Hm are Bessel and Hankel func- 
tions of the first kind, respectively. The resonances /cres 
of the microdisk are given by the poles of the Green's 
function ([1]), i.e. they satisfy the condition 

nHm{kresR)Jm{kresnR) - Hj^{k resR) Jm{k resTlR) = 0. 

Resonances differing by the sign of m are degenerate. 

Microdisks with a point scatterer. — The Green's func- 
tion ([1]) is logarithmically divergent at the point r = tq = 
d where d < R. Since 

Ho{z) = l + i-(lnz + 7-ln2) + 0(z') 

TT 

where 7 = 0.5772156649 ... is the Euler-Mascheroni con- 
stant, the regularized Green's function Gr can be ob- 
tained from ^ if we subtract the term In {ko\r — ro|)/27r, 
where ko is an arbitrary constant. It then follows from 
self- adjoint extension theory Q that the resonances kres 
of the microdisk with a point scatterer at position d and 
of coupling strength A satisfy the condition ^] 

= l-Aa(d,d,A:res), (2) 



where Gr(d,d, /c) is the regularized Green's function at 
the point r = tq = d. It is convenient to intro- 
duce the new coupling parameter a which is defined by 
27r/A = — In {k^a) and hence is a combination of the free 
parameter ko and the coupling strength A. Equation ([2]) 
then reads 

i 1 k TiiQj i C 

The parameter a can be directly interpreted in the limit 
/ca <C 1 as the radius of a localized perturbation of the 

disk [ItI. 

The resonance wavefunctions Ez{t^ ^res) are given by 

Ez (r, kres) =NG{t, d, kres). 

where G(r, d, /cres) is the unregularized Green's func- 
tion ([T|) and A/" is a normalization factor. 

Note that m is no longer a good quantum number of 
the microdisk with a point scatterer. However, it fol- 
lows from the condition ([2]) that for a vanishing strength 
of the point scatterer (A = 0± or equivalently a = 
or a = +00) the resonances of the microdisk with scat- 
terer coincide with those of the unperturbed microdisk, 
which, for m 7^ 0, are twofold degenerate. In fact, if the 
coupling parameter a varies from to +00 one of the 
members of each pair of degenerate resonances remains 
unchanged. The other member is moving in the complex 
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k plane along a line segment that connects two resonances 
of the unperturbed microdisk. This can be understood 
from the fact that for a vanishing strength of the scatterer 
the angular dependence of the resonant modes is given 
by sm{m(p) and cos{m(p). The unperturbed resonance is 
the one with a nodal line along the x-axis on which we 
place the scatterer (see Fig. [T]), i.e. the resonance with 
the angular part s'm{m(p). 

We study the dynamics of the resonances upon varying 
the coupling parameter a for a typical GaAs microdisk of 
effective refractive index n = 3 and radius R = 1 fim [18|, 
with a point scatterer placed at three different distances 
(0.99 /im, 0.495 /im, 0.25 jam) from the center of the disk. 
Figures [2] and [3] show the corresponding resonances in the 
frequency regions u = cRe (k) /27r = 0.048 - 1.098 x 10^"^ 
Hz (mid and near infrared) and u = 5.806 — 6.035 x 
10^^ Hz (green light), respectively. Interestingly, the line 
segments parametrized by the coupling parameter a do 
not only connect different resonances of the unperturbed 
microdisk but as indicated in the insets in Fig. [2] and in 
Fig. [3] there are also loops connecting single resonances 
to themselves. 
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FIG. 2: (color online). Level dynamics of the resonances 
in the complex wavenumber plane for a dielectric disk with 
n = 3.0 and R = 1 /xm and a point scatterer of varying cou- 
pling parameter a. The quantum numbers of the unperturbed 
resonances are marked in brackets. For the upper curve the 
scatterer has distance d = 0.99 /xm from the centre, for the 
middle curve d = 0.495 /xm, and for the lower curve 0.25 /xm. 
The color code indicates the directionality of the emission 
(green marks small values of A|y|2, red marks high values of 
A|y|2). The cross corresponds to the HD mode (A|y|2 = 2.12 
Q 251) with kR 
a ^ 0.754). 



2.0571 - i 0.0164 {d = 0.495 /xm, 



In order to quantify the far-field directionality of the 
electric field we consider its asymptotic behaviour for r 
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FIG. 3: (color online). Continuation of Fig. [2] for a different 
frequency range. The lower cross corresponds to the HD mode 
(A|^|2 = 2.36, Q = 212) with kR = 12.2257 - i0.1156 {d = 
0.25 /im, a ^ 0.047). The upper cross corresponds to the HD 
mode (A|y|2 = 4.71, Q = 15700) with kR = 12. 5513-i 0.0016 
((i = 0.495 /im, a ^ 0.002). 




FIG. 4: Polar plot of the far-field intensity \f {(p)\^ for the 
unperturbed resonant mode (0, 1) (dashed circle), the unper- 
turbed resonant mode (2,2) (thin solid line), and the HD 
perturbed resonant mode marked by the cross in Fig. [2] (thick 
solid line). The parameters are the same as in Fig. [2] 



oo which has the form 

Ez (r, kres) = (r, (p, kres) OC ^^P(^^^es^) ^ ^ 

To characterize the directionality we compute the nor- 
malized variance of the far-field intensity 
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From this definition it follows that A 



and A 



0.5 for resonances which have m = and m 7^ 0, respec- 
tively, for a = or a = 00. 

In Figs. [2] and [3] the directionality A|j|2 is indicated 
by color in the range from light green (low directional 
modes) through blue and violet to deep red (highly di- 
rectional modes). It reaches values as high as A|j|2 ^ 5 
for some specific modes. We see that there are highly 
directional modes (HD) for a wide range of Q-factors 
which in terms of the complex wavenumbers k are given 
by Q = -Re(A:)/(2Im(A:)). 

To illustrate the directionality in more detail we com- 
pare in Fig. m the function \f {(f)f for two unperturbed 





FIG. 5: (color online). The intensity of the electric field in 
near- and far-field regions for unperturbed resonant modes 
(0,1) (left), (2,2) (middle), and for the HD perturbed reso- 
nant mode marked by the cross in Fig. [2] (right). The param- 
eters are the same as in Fig. (2] 



FIG. 7: (color online). The intensity of the electric field in 
near- and far-field regions for HD perturbed resonant mode 
shown by the lower cross in Fig. [3] (left) and for the HD per- 
turbed resonant mode marked by the upper cross in Fig. [3] 
(right). The parameters are the same as in Fig. [3] 




FIG. 6: Polar plot of the far- field intensity | /((/?) |^ for the un- 
perturbed resonance with quantum numbers (m, q) = (5,10) 
(thin solid line in the top panel), the HD perturbed resonant 
mode marked by the lower cross in Fig. [3] (thick solid line 
in the top panel), the unperturbed resonance with {m,q) = 
(21, 4) (thin solid line in the bottom panel), the HD perturbed 
resonant mode marked by the upper cross in Fig. E] (thick solid 
line in the bottom panel). The parameters are the same as in 
Fig. El 



resonant modes and one perturbed HD mode in the fre- 
quency region of Fig. O The corresponding near- and 
far-field electric field intensities are shown in Fig. [5l For 
a vanishing coupling strength the angular dependence is 
sinusoidally modulated with the number of minima be- 
ing given by the quantum number m, i.e. only for m = 
is the emission isotropic. The radial quantum number q 
determines the number of minima in the radial direction. 

Figures [6] and [3 show the near- and far-field electric 
field intensities for two perturbed HD modes in the higher 
frequency region of Fig. [3l This demonstrates that these 
can differ significantly from the nearby modes of the un- 
perturbed system also shown in Fig. [H While the near- 
field pattern of the mode shown in the right panel of 
Fig. [71 is still reminiscent of a WGM, the mode in the left 
panel is strongly scarred. The far-field patterns shown 
in Figs, ini and (3 also show that the scatterer can cause 
directional output in either direction of the symmetry 
axis. The mode in the bottom panel of Fig. |6l has both 
extremely high directionality and a very high Q-factor. 
From Fig. |3l we see that there is quite a broad spectral 
range of green light with these properties the realization 
of which is a major goal in semiconductor physics. 

Conclusions. — In summary, we presented theoretical 
results that demonstrate the existence of highly direc- 



tional TM-modes in the emission spectrum of a microdisk 
cavity with a point scatterer. These modes can appear 
even for a scatterer with a very weak coupling constant 
which promises the feasibility of an experimental real- 
ization of such cavities. The level dynamics of the reso- 
nances upon varying the coupling constant is of theoret- 
ical interest in its own right and deserves further investi- 
gation. Similarly, it would be interesting to get a deeper 
insight into the output directionality by relating the res- 
onance wavefunctions to the underlying ray dynamics in 
the semiclassical limit. 
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